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Abstract 

We applied canonical transformation to water wave equation not only to remove 
cubic nonlinear terms but to simplify drastically fourth order terms in Hamiltonian. 
This transformation explicitly uses the fact of vanishing exact four waves interaction 
for water gravity waves for 2D potential fluid. After the transformation well-known 
but cumbersome Zakharov equation is drastically simplified and can be written in 
X-space in compact way. This new equation is very suitable as for analytic study 
as for numerical simulation. Localized in space breather-type solution was found. 
Numerical simulation of collision of two such breathers strongly supports hypothesis 
of integrability of 2-D free surface hydrodynamics. 
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1 Introduction 



The work described here is motivated by two remarkable facts regarding one- 
dimensional free surface hydrodynamics: 

• In [lj it was shown that four-wave interaction coefficient vanishes on the 
resonant manifold 
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k + ki = k 2 + k 3 , 

• In [2113] it was demonstrated that giant breather, highly nonlinear, exists on 
the fluid surface without radiation. Moreover, space-time spectrum of the 
breather consists of waves propagating in the same direction. 

These two facts bring the idea that fourth order wave interactions can be dras- 
tically simplified by some canonical transformation of the Hamiltonian. Below 
we show how this transformation looks like. Dynamic equation derived using 
this transformation is very elegant and simple. It can be easily generalized for 
"almost" one- dimensional waves. 



2 Compact equation 



A one- dimensional potential flow of an ideal incompressible fluid with a free 
surface in a gravity field fluid is described by the following set of equations: 



4>xx + 4>zz = 
Vt + Vx4>x = 0; 



— > 0, Z — > — oo), 
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here rj(x, t) - is the shape of a surface, <j)(x, z, t) - is a potential function of the 
flow and g - is a gravitational acceleration. As was shown in|3], the variables 



f](x,t) and ip(x,t) = <p(x,z,t) 
equations 



are canonically conjugated, and satisfy the 



z=r] 



dip 

~dt 



6H_ 

5rj 



Or] 

at 



SH_ 



Here H = K + U is the total energy of the fluid with the following kinetic and 
potential energy terms: 



K — — I dx I v 2 dz U = ^ / r] 2 dx 
It is convenient to introduce normal complex variable a^: 
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Vk = J^(a k + a*_ k ) ^ k = -i^^-(a k - a*^ k ) 

here u k = y/gk -is the dispersion law for the gravity waves, and Fourier trans- 
formations ip{x) — > ipk and r/(x) — > r] k are defined as follows: 

fk = -4= / f(x)e~ tkx dx, f{x) = ^=1 f k e +ik »dk. 

Hamiltonian can be expanded in an infinite series in powers of (see [lf5] ) 



H = H 2 + H 3 + HA+ ... 
This variable satisfies the equation 



<9a fc 8H 
at oa* k 



where 



H 2 = J u k a k a* k dk, 

H 3 = J V k k ik2 {a* k a kl a k2 + a k a* ki a* k2 }5 k - kl ^ k jlkdk 1 dk2 
+ ^ J U kklk2 {a k a kl a k2 + a* k a* ki a* k2 }5 k+kl+k . 2 dkdk 1 dk 2 
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L kkl = (kki) + |fc||fci| 
Fourth order part of Hamiltonian is the following: 



H 4 = - J W k ^ k2 al i al 2 a k:j a k4 5 kl+k2 ^ k:i ^ k4 dk 1 dk 2 dk 3 dk^+ 

+ g J G k klk2ks {a ki a k2 a* hi a k4 + a kl a k2 a k3 a* k4 5 kl+k2+k:i ^ k4 dkidk 2 dk 3 dk 4 + 

+ — J Rk 1 k2k a k4( a k 1 a k2 a k: i a k4 + a k i ak 2 a k s ak i )bk 1 +k 2 +k: i +k i dkidk 2 dk 3 dk4 : 
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Here W^, Gfr lkah3 and R kl k 2 k 3 k 4 are equal to: 
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Here 

M S 2 4 = |fc x A; 2 | 5 |A; 3 fc 4 |i(|A;i + A;3| + |A;i + A;4| + |fc2 + /c3| + 1^2 + ^1 -2|A;i|-2|A;2|). 

Now one can apply canonical transformation from variables a/c to b k to exclude 
non resonant cubic terms along with non resonant fourth order terms with 
coefficients G k ^ ks and Rk 1 k 2 k 3 k 4 - This transformation up to the accuracy 0(b 5 ) 
has the form [4"f6fT] : 

d'k = bk + J F^ ik2 bk 1 bk 2 5k-k 1 -k 2 dkidk2 — 2 / T k z ki bl i bk 2 5k+k 1 -k 2 dk 1 dk 2 + 
+ J T kklk2 b* ki b* k2 5 k+kl+k2 dk 1 dk 2 + J B k 2 k k:i b* ki b k2 b k . A 5 k +k 1 -k 2 -k 3 dk 1 dk 2 dk 3 + 
+ J C k l lk2 bl 1 bl 2 b ka 5 k+kl+k2 - k3 dk 1 dk 2 dk 3 + J Sk^k^bl^bljk+k^+^dhdkidkl) 

Here 
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and B k k ^ is an arbitrary function satisfying the following symmetry condi- 
tions: 



TDkki Dfcifc TDkki ( TDk 2 k 3 \* 

n k 2 k 3 — n k 2 k 3 — n k 3 k 2 — \ D kkx ) ■ 

Coefficients C kk k and S k kik 2 k 3 provide vanishing corresponding forth-order 
terms in the new Hamiltonian. 

Details of this transformation can be found in [515] . In K-space Hamiltonian 
has the form: 
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U = J u k \b k \ 2 + - J f^bl bkzbkabkiSfa+kz-ks-ktdkidkzdksdkt 
In X-space it corresponds to: 



n = J b*u k bdx + 



,d_ 

dx 



2 d 



d, - ^ \ \b\* ■ k{\bfim) 



Here b' = Jj^b. After integrating by parts Hamiltonian acquires very nice form: 



H = Jb*u k bdx + ^J\b'\ 2 



(W* - 6*6') - K\b\' 



dx. 



(2.3) 



Corresponding equation of motion is the following: 



db i 

l m =Ukb+ 8 



b*^-(b' 2 ) - #(6*'^ 



dx 



dx dx 



b-K(\bf)-^(b'k(\bm 
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3 Monochromatic wave and modulational instability 

Monochromatic wave 

b{x) = B e i{ - kox -"^ (3.5) 

is the simplest solution of (12 ,4p . Indeed, plugging (13. 5 p in to the equation (12. 4p 
one can get the following relation 

uq = u) ko + -kl\B Q \ 2 . (3.6) 

Recalling transformation from a k to b k one can see that for waves with small 
amplitude ( a k ~ b k ) 

\Bo\ -j^Vo, 

and relation (I3.6P coincides with well known Stokes correction to the frequency 
due to finite wave amplitude. 

Uo = u ko (l + t^M 2 )- ( 3 - 7 ) 
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B = -L / h^-^dx 
V In J 



3. 1 Modulational instability of monochromatic wave 

Let us consider perturbation to the solution 

b = B e i{koX -" ot) 

where 

B n = - 

/2?r 

and 

= u; fc0 + \\B,\ 2 kl ^\b k0 \ 2 kl = \T k £\b k0 \ 2 . 
Perturbed solution has the following form: 

b => (b k0 + 5b k0+k e- inkt + Sh^e-^e-^. (3. 

with the following condition: 

Vl k = — 

Plugging perturbed solution ( 13. 8 p in to the equation 

1 f ~ 

ib k = u k b k + - J T^ k \bl 1 bk 2 bk 3 5k+k 1 -k 2 ~kjlk 1 dk'rfk3 
we get the sum of two independent equations: 

i8b ko+k + (u + Q k )8b ko+k = u ko+k 5b ko+k + T k °+ kk °\b ko \ 2 5b ko+k + -f k ° k + o kk ° k b 2 kQ 5b* ko _, 
+ i8bk -k + (w + Q- k )5b ko ^ k = u ko - k 6b ko - k + f k °I kk °\b ko \ 2 5b ko ^ k + - 

Expressions for f k °l kk ° and T k ° k kko+k can be easily obtained: 



isi^-k + (u + n_ k )5b k0 ^ k = u k0 „ k 5b k ^ k + f k °: kk °\b k0 \ 2 5b k0 „ k + lf k ° k0 kk0+k b 2 k M 0+ k 



e — iu> t— if2fet_|_ 



e —iu>ot-iU_ k t 



= | + H \ n * + ^ - ko \k\ + k % (3.9) 

r i"o fc/£o+/£ = ^o 2 -^o|^|-y). (3.10) 
Looking at even and odd powers of k one can see that 

n k = u * +h - u *- h + ^(u*-\k\)k 

Let us denote 



Then 

iSb ko+k = d(k)5b ko+k - 



2 



Suppose Sb ko+k growth as 

Sb ko+ k 5b ko+k e lkt 
one can easily obtain the following formula for ^ k : 

31 S 



ll 



-d(k) - 



-k k 2 



d(k) + \B \ 2 k (k - 1 -^-) 2 



(3.11) 



If we introduce steepness of the carrier wave u ko fi = \B \ k^ and approximate 
d{k) as 



d{k) 



1 1 k 2 

— » = Wk-r, 

8 fc ° 8 °& 2 



then for growth rate is 



Ioj 



-^(l-6/x 2 )fc 2 
/c 



H 2 (,k 



|k| 

T 



(3.12) 



The difference between this formula and well-known expression derived from 
the nonlinear Schrodinger equation is highlighted by two boldfaced terms. 



4 Breathers and numerical sumulation of its collisions 

Breather is the localized solution of (12. 4p of the following type: 

b(x,t) = B(x - vt)e i{kox -" ot \ (4.13) 

where ko is the wavenumber of the carrier wave, V is the group velocity and 
Co>o is the frequency close to u ko . In the Fourier space breather can be written 
as follow: 

h(t) = e-^ m+Vk ^ k , (4.14) 

where f2 is close to ^p-. 

For (fi k the following equation is valid: 

(Q + Vk - u k )(f) k = J f kkuk2k3 (pl i (j) k2 (p k .J k+kl ^ k , 2 ^ k , A dkidk2dk 3 . (4.15) 

One can treat <p k as pure real function of k. 
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Fig. 1. Real part of b(x) with V = 0.1 and U = 2.53. 
To solve equation (14. 15j) one can use Petviashvili iteration method: 



(ft + V k - w fc )#b = Mn J nk^Jty^lSk+^^dhdhdh. (4.16) 
Petviashvili coefficient M is the following: 

<<ft(Sl + Vk-u> k )ffi> 



_, 3 

2 



Below we present typical numerical solution of (I4.15p . Calculation were made 
in the periodic domain 2ir with carier wavenumber k Q ~ 25, V = 0.1 and 
= 2.53EH In the Figures [U |2l |3] one can see real part offr(x), modulus of 
b(x) and Fourier spectrum of b(x). 

Very important question from the point of view of integrability of the equa- 
tion (12.41) is the question about collision of two breathers. To study breathers 
collision we performed the following numerical simulation: 

• As initial condition we have used two beathers separated in space (distance 
was 7T.) 



1 steepness of carrier wave fi must not exeed the crical value y | as it follows from 
formula for growt rate of modulatinal instability (|3.12p 
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Fig. 2. Modulus of b(x) with V = 0.1 and VL = 2.53. (Recall envelope in the Nonlinear 
Shrodinger Equation 
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Fig. 3. Spectrum of b(x) with V = 0.1 and O = 2.53. 
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Fig. 4. Initial condition with two breathers. 

5.1, Vi 



0.05. Carrier 



• First breather has the following parameters: fii 
wave number appears to be ~ 100. 

• For the second breather - Q 2 = 2.53, V2 = 0.1. Carrier wave number appears 
to be ~ 25 

This initial condition is show in Figure HI Its Fourier spectrum is shown in 
Figure After time Ty^o — 62.8 breathers collide. In the Figures E] and [7] 
one can see breathers at the time close to collision (t = 50) and at the moment 
of collision (t = 63). Fourier spectrum of two breathers at t = 63 is shown 
in Figure [H And finally we show the picture of two breathers at t — 126 
when they separated again at distance ~ it. Real part of b(x) and Fourier 
spectrum of that is given in Figures [91 [T0J So, the simulation demonstrates 
no interaction of breathers. It suggests that equation (12 .4p is integrable. 



5 Conclusion 



Simple equation describing evolution of 1-D water waves is derived. Derivation 
of this equation is based on the important property of vanishing four-wave in- 
teraction for gravity water waves. This property allows to simplify drastically 
well-known Zakharov's equation for water waves, which is very cumbersome. 
Written in X-space instead of K-space, it allows further analytical and numer- 
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Fig. 5. Fourier spectrum of two breathers. 
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Fig. 6. Two breathers just before collision. 
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Fig. 7. Two breathers collides. 
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Fig. 8. Fourier spectrum at the moment of collision. 
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Fig. 9. Two breathers after collision. 
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Fig. 10. Fourier spectrum at the moment of collision. 
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ical study. The equation has breather-type solution which was found numeri- 
cally using Petviashvili iteration method. Numerical simulation of collision of 
two breathers shows behavior which is typical for integrable system. It can be 
considered as numerical proof of integrability. 

This new equation can be generalized for the "almost" 2-D waves, or "almost" 
3-D fluid. When considering waves slightly inhomogeneous in transverse di- 
rection, one can think in the spirit of Kadomtsev-Petviashvili equation for 
Korteveg-de-Vries equation, namely one can treat now frequency u k as two 
dimensional, tOk x ,k y , while leaving coefficient not dependent on y. b now 
depends on both x and y: 



n = J b*u kx , ky bdxdy +\J \b' x \ 2 [i(M£ - b%) 



K x \b\" 



dxdy. (5.17) 



6 Acknowledgments 



This work was supported by Grant of Government of Russian Federation for 
support of scientific research, carried under direction of leading scientists in 
russian universities Nil. G34. 31. 0035 (leading scientist - Zakharov V.E., GOU 
VPO "Novosibirsk State University"). Also was it was supported by the US 
Army Corps of Engineers Grant W912-BU-08-P-0143, by ONR Grant N00014- 
10-1-0991, NSF Grant DMS 0404577, Grant NOPP "TSA-a two scale approx- 
imation for wind-generated ocean surface waves", RFBR Grant 09-01-00631 
and RFBR Grant 09-05-13605, the Program "Fundamental Problems in Non- 
linear Dynamics" from the RAS Presidium, and Grant "Leading Scientific 
Schools of Russia" . 



References 



[1] Dyachenko A.I. and V.E.Zakharov V.E., Phys. Lett. A 190, (1994), 144-148. 

[2] Dyachenko A.I. and Zakharov V.E, Pis'ma v ZhETF, 88, (2008) 356. 

[3] Zakharov V.E and Dyachenko A.I, European Journal of Mechanics, B/Fluids, 
29(2), (2010), 127-131. 

[4] Zakharov, V.E. , Prikl.Mekh.Tekh.Fiz.(in Russian), 2, (1968), 190 

[5] Crawford,D.E., Yuen,H.G. and Saffman,P.G., Wave Motion, 2, (1980), 1 

[6] Krasitskii,V.P., Sov. Phys. JETP, 71,(1990), 921 



14 



[7] Zakharov V.E., Lvov V.S. and Falkovich G., "Kolmogorov Spectra of 
Turbulence I", Springer- Verlag, 1992. 

[8] A.I. Dyachenko and V.E. Zakharov, "Compact equation for gravity waves on 
deep water", JETP Letters, 93, (12), (2011) 701. 

[9] A.I. Dyachenko, V.E. Zakharov, "On dynamical equation for water waves in one 
horizontal dimension", European Journal of Mechanics - B/Fluids, In Press, 
Accepted Manuscript, Available online 7 September 2011. 



15 



